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Abstract 

In  this  article,  we  discuss  techniques  for  sense  through  wall  human  detection  for  different  types  of  walls.  We  have 
focused  on  detection  of  stationary  human  target  behind  wall  based  on  breathing  movements.  In  detecting  the 
breathing  motion,  a  Doppler  based  method  is  used.  Also  a  new  approach  based  on  short  time  Fourier  transform  is 
discussed  and  an  already  proposed  clutter  reduction  technique  based  on  singular  value  decomposition  is  applied 
to  different  measurements. 
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I  Introduction 

Detection  of  human  target  through  wall  is  of  interest  for 
many  applications.  Military  industry  could  use  it  for 
hostage  rescue  situations.  In  such  scenarios,  detection 
and  location  of  humans  inside  a  room  is  very  critical  as 
unknown  building  layout  together  with  presence  of 
armed  persons  can  be  dangerous  for  the  rescuers. 
Another  use  could  be  for  disaster  search  and  rescue 
operations  such  as  people  trapped  under  building  debris 
during  earthquake,  explosion  or  fire. 

Ultra  WideBand  (UWB)  technology  has  emerged  as 
one  of  the  preferred  choices  for  such  applications  due 
to  its  good  range  resolution  and  good  penetration 
through  most  of  the  building  materials.  High  range  reso¬ 
lution  is  a  result  of  high  bandwidth  of  UWB  radar  and  it 
helps  in  better  separation  of  multiple  targets.  Detection 
of  human  target  is  based  on  the  fact  that  there  is  always 
some  movement  due  to  breathing  or  movement  of  body 
parts  (as  in  case  of  a  walking  person).  This  small  move¬ 
ment  can  be  used  to  detect  a  human  being  from  other 
objects  behind  a  wall  or  beneath  rubble  but  it  becomes 
challenging  due  to  high  clutter  from  the  wall  and  other 
objects  inside  a  room. 
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The  focus  of  this  article  is  on  detection  techniques  for 
a  motionless  human  target  using  a  monostatic  UWB 
radar. 

II  UWB  Overview  and  Features 

UWB  systems  are  the  ones  which  use  signals  with  mini¬ 
mum  (10  dB)  bandwidth  of  500  MHz  or  fractional 
bandwidth  of  at  least  20%. 

-  fT 

Fractional  bandwidth  =  2  - - —  (1) 

/h+/l 

where, /h  and  /L  are  highest  and  lowest  frequency 
points,  respectively,  with  signal  10  dB  below  peak 
emission. 

A.  Large  bandwidth-high  range  resoiution 

The  relation  between  pulse  width  and  radar  range  reso¬ 
lution  is  given  as 

T  •  C  C 

Range  resolution  =  — =  —  (2) 

where  x  =  Pulse  width  in  time  domain,  B  =  bandwidth 
of  the  pulse,  and 

c  =  speed  of  electromagnetic  waves. 

Good  range  resolution  property  of  UWB  can  be  used 
for  localizing  the  target  in  an  indoor  environment. 
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B.  UWB  radar  penetration  through  wall  [1] 

As  per  the  Electromagnetic  theory,  lower  frequencies 
have  better  penetrating  properties.  UWB  radar  uses  a 
large  spectrum  in  combination  with  lower  frequencies 
which  makes  it  suitable  for  applications  such  as  ground 
penetrating  radar,  foliage  penetrating  radar  [2],  and 
short-range  radar  to  detect  hidden  objects  behind  walls. 
This  penetration  property  is  also  of  great  importance  for 
indoor  location  systems.  Shorter  wavelength  makes  pos¬ 
sible  use  of  smaller  dimensions  of  receive  and  transmit 
antennas.  On  the  other  hand,  an  increase  in  a  center 
wavelength  of  the  signal  is  desirable  for  enhancing  the 
penetrating  capability  of  electromagnetic  waves  through 
walls.  However,  an  increase  in  the  wavelength  is  again 
restricted  by  two  factors:  the  first  one  is  related  with 
shielding  sounding  signals  by  metallic  meshes  in  con¬ 
crete  walls,  while  the  second  one  decreases  the  RCS  of 
the  target  when  the  wavelength  exceeds  the  sizes  of  the 
target.  The  estimation  carried  out  have  shown  that  for 
conducting  of  rescue  activities  in  ruins,  typical  concrete 
buildings  and  facilities  the  most  optimal  is  the  frequency 
range  0.8  to  2  GHz. 

Ill  Effects  of  wall  and  human  body  as  radar  target 

A.  Wall  clutter  [3,4] 

In  through  wall  target  detection,  clutter  can  be  due  to 
many  reasons  like  wall  coupling,  antenna  coupling,  mul¬ 
tiple  reflections. 

In  through  wall  target  detection,  clutter  reduction  can 
play  important  part  to  accurately  detect  the  target  and 
remove  the  unwanted  signals  which  arise  due  to  the 
reflection  from  the  wall  and  other  reflections  due  to 
unwanted  objects.  Once  the  signal  is  transmitted 
through  the  antenna,  it  suffers  attenuation  due  to  wall 
and  other  obstacles.  A  clutter  reduction  technique  such 
as  SVD  reduces  signal  due  to  wall  and  enhances  the 
peak  due  to  target. 

B.  Human  target  detection  [4-6] 

Detection  of  human  beings  with  radars  is  based  on 
movement  detection  (e.g.,  walking  human),  chest  move¬ 
ments  due  to  breathing  or  heartbeat.  Heart  beat  and 
respiratory  motions  cause  changes  in  frequency,  phase, 
amplitude,  and  arrival  time  of  reflected  signal  from  a 
human  being.  In  case  of  through  wall  human  target 
detection,  these  changes  can  be  very  small,  especially  for 
a  brick  or  concrete  walls. 

Reflected  UWB  signal  is  highly  sensitive  to  human 
posture  and  thus  makes  detection  process  challenging. 
For  example,  the  signal  reflected  from  the  breathing 
human  causes  changes  in  received  waveform  shape. 

An  effective  human  detection  method  requires  a 
model  of  UWB  radar  waveform  propagation  and  scatter¬ 
ing,  e.g.,  interaction  with  the  human  body.  A  perfectly 


reflecting  target  e.g.  a  metal  plate  with  an  infinite  area 
returns  the  incident  UWB  pulse  along  a  single-path. 
However,  for  a  target  such  as  human  body,  which  has 
complex  shape  and  whose  spatial  extent  is  larger  than 
the  transmitted  UWB  signal  pulse  width,  the  returned 
UWB  radar  signal  consists  of  multipath  components,  as 
the  incident  UWB  pulse  scatters  independently  from  dif¬ 
ferent  human  body  parts  at  different  times  with  different 
amplitudes  (depending  on  the  distance  to  the  body  part 
and  the  size,  shape,  and  composition  of  the  scattering 
part). 

IV.  Measurements 

This  article  considers  P220  UWB  radar  in  monostatic 
mode  (shown  in  Figure  1)  where  waveform  pulses  are 
transmitted  from  a  single  Omni-directional  antenna  and 
the  scattered  waveforms  are  received  by  a  collocated 
Omni-directional  antenna  [1,7].  The  two  antenna  ports 
on  the  P220  are  used  for  the  transmit  and  receive 
antennas.  An  Ethernet  cable  is  used  to  connect  the 
radio  to  the  PC  and  radar  can  be  controlled  using  appli¬ 
cation  software  provided  with  the  radios.  The  P220 
UWB  radar  used  here  has  center  frequency  of  4.3  GHz 
with  a  10-dB  bandwidth  of  2.3  GHz.  This  radar  provides 
a  resolution  of  6.5  cm. 

In  this  section,  we  look  at  few  important  related  para¬ 
meters  related  to  radio  configuration.  These  parameters 
are  important  in  analyzing  captured  scans. 

Integration  is  the  number  of  radio  pulses  that  radar 
combines  to  increase  the  signal-to-noise  ratio.  It  is  the 
total  number  of  UWB  pulses  per  waveform  (scan) 
sample. 

Window  Size  (ft)  is  the  width  of  the  window’,  in  which 
motion  can  be  detected. 

Pulses  Per  Waveform  is  the  number  of  UWB  radio 
pulses  required  for  the  entire  waveform  (single  scan). 
Divide  this  by  the  pulse  rate  to  determine  the  theoretical 
maximum  scan  rate. 

Step  Size  (ps)  is  the  waveform  scan  resolution  (step 
size  between  points),  in  picoseconds  (1  bin  =  3.18  pS). 
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A.  Data  collection 

For  each  measurement  set,  scans  were  acquired  for 
duration  of  around  1  min.  The  number  of  scans 
acquired  depends  on  the  scan  rate  which  in  turn 
depends  on  the  waveform  scan  resolution,  the  window 
size,  and  the  Integration  size.  Also  scans  were  taken 
once  with  human  target  and  once  without  human  target. 

The  equations  below  show  the  relation  between  some 
important  parameters. 

.  _  2  *  window  size  *  0.3049 

Number  or  sample  points  per  scan  = - - - - -  (J) 

c  *  step  size  (bins) 

Number  of  pulses  per  scan  =  Integration* number  of  sample  points  per  scan  (4) 

Number  of  pulses/scan 

Scan  rate  =  — - — - 7 -  (5) 

Pulse  repetition  frequency 

Total  number  of  scans  collected  =  scan  rate  *  total  data  collection  time  (6) 

From  the  above  expressions,  we  can  see  that  increas¬ 
ing  the  scan  Window  Size  or  Integration  size  increases 
the  scan  time  and  thus  reduces  the  Scan  Rate .  However, 
increasing  the  Step  Size  decreases  the  Scan  Rate . 

B.  Measurement  locations 

For  the  purpose  of  this  project,  measurements  were  taken 
at  four  different  locations  having  different  types  of  walls. 
The  radar  parameters  in  each  of  the  cases  given  below  were 
Integration:  Hardware  Integration  =  512,  Software 
Integration  =  2,  Pulse  Repetition  Frequency:  9.6  MHz 
Step  Size:  1  bin,  7  bin,  Window  Size  (ft):  10  ft 

(1)  Gypsum  wall 

Figure  2  shows  the  location  of  the  radar  and  Human 
target  on  different  sides  of  a  1-ft  thick  Gypsum  partition 
wall.  Person  is  at  a  distance  of  6.5  ft  from  the  radar  on 
the  other  side  of  the  wall  and  the  height  of  the  antennas 
from  ground  is  3’4”. 


(2)  Wooden  door 

Figure  3  shows  the  location  of  the  radar  and  Human 
target  on  different  sides  of  a  4-cm  wooden  door.  Person 
is  standing  at  a  distance  of  7’6”  from  the  radar  on  the 
other  side  of  the  door  and  the  height  of  the  antennas 
from  ground  is  3’4  “. 

(3)  Brick  wall 

Figure  4  shows  the  location  of  the  radar  and  Human 
target  on  different  sides  of  a  12-cm  Brick  wall.  Person  is 
standing  at  a  distance  of  8’  from  the  radar  on  the  other 
side  of  the  door  and  the  height  of  the  antennas  from 
ground  is  3’4”. 

(4)  Load  bearing  concrete  wall 

In  this  case,  measurements  were  taken  at  two  different 
positions  as  shown  in  Figure  5.  In  both  cases,  person  is 
standing  at  7’6”  from  the  radar  and  the  height  of  the 
radar  is  3’4”. 

V.  Measurement  analysis 

In  this  section,  we  discuss  the  three  approaches  that  are 
used  in  this  article  to  analyze  the  measurements. 

A.  Detection  of  breathing  movements 

This  approach  is  based  on  detection  of  small  chest 
movements  associated  with  a  breathing  motionless 
human.  This  motion  is  very  small  and  results  in  very 
weak  radar  echo.  However,  since  it  is  periodic  motion 
it  can  be  detected  by  application  of  signal  processing 
techniques  which  enhances  the  ‘breathing’  signal  from 
noise. 

Breathing  motion  will  cause  periodic  changes  in  the 
received  signal  at  a  distance  where  target  is  located. 
This  periodic  change  is  reflected  across  multiple  scans. 
Thus  an  N  x  M  matrix  A  is  constructed  using  ‘M’ 
scans,  each  of  length  ‘AT,  as  columns  of  matrix  A.  Then 
difference  is  taken  between  successive  columns  of 
matrix  A,  which  captures  changes  from  one  scan  to 
another  and  helps  to  suppress  the  static  clutter  signal. 
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Finally  DFT  is  performed  on  each  row  of  the  resulting 
matrix  which  clearly  shows  the  breathing  human  target. 
This  approach  is  summarized  below. 

Step  1.  Matrix  A  constructed  using  ‘AT  scans  arranged 
in  columns 

~  Scan  1  Scan  2  Scan  3  •  •  •  Scan  Af 

sample  1  sample  1  sample  1  •  •  •  sample  1 
^  _  sample2  sample2  sample2  •  •  •  sample2  ^ 

_  sample  N  sample  N  sample  N  •  •  •  sample  N  _ 

Step  2.  Matrix  D  is  the  difference  between  successive 
columns  of  A 


"Scanl  -  Scan2 

Scan2  -  Scan3 

Scan(A4  -  1)  -  ScanA/T 

samplel 

samplel 

samplel 

D  = 

sample2 

sample2 

sample2 

(8) 

_  sampleN 

sampleN 

sampleN 

Step  3.  Take  Discrete  Fourier  Transform  of  each  row 
of  the  Matrix  D. 

This  technique  works  for  gypsum  wall,  wooden  door, 
and  brick  wall.  Below  are  the  observations  for  these 
cases. 


(1)  Gypsum  wall 

Figure  6  shows  the  D-Matrix  with  and  without  target.  In 
this  case,  D  matrix  is  constructed  using  100  scans  cap¬ 
tured  at  scan  rate  of  0.6827  scans/s  for  total  time  dura¬ 
tion  of  68  s. 

A  discrete  Fourier  transform  (DFT)  on  each  row  of  D- 
matrix  shows  the  breathing  rate  of  a  human  target  at  6.5 
ft.  (Figure  7). 

Figure  8  shows  the  case  where  the  person  is  moving 
his  hands  towards  the  radar  and  back  at  rate  close  to  1 
Hz 

(2)  Wooden  door 

See  Figures  9  and  10. 

(3)  Brick  wall 
See  Figure  11. 

This  method  of  detecting  motionless  people  may  not 
work  in  all  cases.  For  example,  this  approach  works  well 
for  wooden  door,  gypsum  partition  wall,  and  brick  wall 
as  shown  above  but  fails  when  the  attenuation  for  signal 
scattered  from  target  is  large  compared  to  the  signal 
reflected  from  wall  or  other  stationary  objects  (e.g.,  con¬ 
crete  wall).  In  such  cases,  detection  of  weak  target  signal 
in  presence  of  strong  clutter  from  wall  is  difficult  and 
will  require  use  of  some  kind  of  clutter  reduction 
method.  Also  this  method  may  fail  when  the  person  has 
his  back  towards  the  wall  as  the  chest  movements  may 
not  be  captured  in  the  resulting  scans. 
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B.  Clutter  reduction  using  SVD  [8,9] 

SVD  is  used  here  to  reduce  wall  clutter.  The  main  aim 
of  SVD  is  to  split  the  Scan-Matrix  into  subspaces  which 
correspond  to  clutter,  target,  and  noise  so  that  the  clut¬ 
ter  can  then  be  rejected.  The  Scan-Matrix  is  constructed 
by  arranging  lM  scans  each  of  length  ‘AT  in  matrix  for¬ 
mat  giving  an  N  x  M  Matrix  A.  Each  column  of  this 
matrix  is  a  single  scan  of  length  M.  The  SVD  of  A  is 
given  as 

A  =  USVT  (9) 

where  UTU  =  I;  VTV  =  I;  the  columns  of  U  are  ortho¬ 
normal  eigenvectors  of  AAT,  the  columns  of  V  are 
orthonormal  eigenvectors  of  ATA,  and  S  is  a  diagonal 
matrix  containing  the  square  roots  of  eigenvalues  from 


U  or  V  in  descending  order. 

A  =  OiU\v\  +  cr2U2t,2  +  a3u3 vj  +  . . .  (10) 

A  =  Ml  +M2  +  M3  +  ...  (11) 

where,  M*  is  called  as  the  it h  Eigen-image  of  A 
It  has  been  found  experimentally  that  first  Eigen- 
image  corresponds  to  clutter,  second  Eigen-image  corre¬ 
sponds  to  target  and  the  rest  are  noise  [2].  Therefore, 
we  have 

A  =  Mciutter  +  Vltarget  +  MnoiSe  (12) 

where,  Mciutter  =  <7\U\v\,  Mtarget  =  crzUjV^,  and  Mnoise  =  03W3V3  +  04M4  +  ...  (13) 
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Figure  6  A-matrix  without  target(top)  and  with  target(bottom) 
for  gypsum  wall. 
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This  technique  does  not  work  for  the  case  of  concrete 
wall  and  wooden  door. 

(1)  Gypsum  wall 

Here  after  applying  SVD  to  A-Matrix,  clutter  is  reduced 
and  target  can  be  detected.  Figures  show  the  A-matrix 
with  target  (Figure  12),  Eigen  Image  corresponding  to 
clutter  (Figure  13)  and  target  (Figure  14). 

(2)  Brick  wall 

See  Figures  15,  16  and  17. 


C.  Short  time  Fourier  transform  and  singular  value 
decomposition 

Short  time  Fourier  transform  (STFT)  is  a  tool  to  analyze 
frequency  contents  of  signals  that  vary  in  time.  STFT 
maps  a  signal  into  a  two-dimensional  function  of  time 
and  frequency.  It  represents  a  kind  of  compromised 
view  of  signal  in  time  and  frequency.  However,  this 
information  is  obtained  with  limited  precision  and  this 
precision  is  determined  by  the  window  size. 

Analysis  using  STFT  involves  choosing  appropriate 
window  size  so  as  to  get  good  resolution  in  both  time 
and  frequency  domain  as  there  is  always  a  trade-off 
between  the  two.  The  window  type  is  selected  according 
to  compromise  between  sidelobe  attenuation 
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requirements  and  mainlobe  width  (usually  Hanning  or 
Blackman  are  used).  Window  size  needs  to  be  adapted 
to  the  signal  and  the  information  one  is  looking  for. 

STFT  of  a  single  scan  will  provide  information  content 
about  frequencies  across  the  scan  duration.  Then  SVD  is 
done  on  the  STFT  output  to  see  if  the  target  can  be 
identified  based  on  its  frequency  content. 

However,  selecting  an  STFT  window  size  that  will 
result  in  enough  resolution  in  both  time  and  frequency 
to  identify  the  target  and  its  distance  is  quiet 
challenging. 

Various  window  sizes  were  tried  for  STFT  to  see  if 
there  is  any  difference  in  the  singular  values  obtained 
from  the  SVD  for  target  and  no  target  case. 

(1)  Gypsum  wall 

Figure  18  shows  the  STFT  for  the  no  target  scan  and 
Figure  19  shows  STFT  for  target  case.  Window  size 


used  in  this  case  is  512  with  an  overlap  of  128.  Target  is 
located  around  4155  sample  index. 

The  singular  values  obtained  from  the  STFT  data  are 
normalized  by  dividing  each  value  with  the  maximum 
singular  value,  as  plotted  in  Figures  20  and  21.  It  is 
observed  in  that  there  is  relative  increase  in  the  second 
singular  value  in  case  when  target  is  present.  This  rela¬ 
tive  increase  is  around  0.2.  However,  this  is  not  consis¬ 
tent  when  applied  to  other  cases  of  wooden  door,  brick 
wall,  and  concrete  wall. 

VI.  Conclusion  and  future  work 

For  detection  of  human  target  using  UWB  radar,  various 
sets  of  measurements  were  taken  using  monostatic  radar 
mode.  Data  were  collected  for  different  types  of  walls 
and  doors.  The  scans  collected  were  analyzed  using 
three  different  approaches.  It  is  observed  that  the  heart 
beat  detection  using  Doppler  approach  works  for  woo¬ 
den  door,  gypsum,  and  brick  wall  but  fails  in  case  of  a 
thick  concrete  wall.  A  second  method  using  singular 
value  decomposition  was  used  to  reduce  clutter  and  this 
works  for  brick  and  gypsum  wall  but  again  fails  for  con¬ 
crete  wall  case.  Finally,  we  tried  an  STFT  and  SVD 
method  based  on  the  idea  that  the  received  signal  in 
case  of  presence  of  target  will  result  in  difference  in  fre¬ 
quency  response  compared  to  no  target  case.  In  this 
method,  selection  of  window  size  and  overlap  size  is  a 
challenging  task.  By  applying  SVD  to  the  STFT  output 
it  is  observed,  in  case  of  gypsum  wall,  that  the  second 
singular  value  changes  relatively  in  presence  of  target. 


0  10  20  30  40  50  60 

Time  in  sec 


Figure  13  Eigen  image  of  clutter  for  gypsum  wall. 
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Figure  20  Normalized  singular  values  in  absence  of  target  for 
gypsum  wall. 


Figure  21  Normalized  singular  values  in  presence  of  target  for 
gypsum  wall. 
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However,  the  approach  does  not  work  for  other  cases. 
The  breathing  rate  detection  approach  does  work  for  all 
except  the  concrete  wall  case,  but  cannot  be  entirely 
relied  upon  as  the  detection  process  is  based  on  changes 
in  the  received  signal  amplitude.  For  example,  the 
received  signal  amplitude  may  be  very  weak  for  person 
standing  with  his/her  back  to  the  wall  under 
consideration. 

Given  the  drawbacks  of  different  approaches  used  in 
this  project,  a  new  approach  is  required  which  should 
work  for  all  cases.  Also  while  collecting  data,  only 
monostatic  mode  of  operation  was  considered.  There¬ 
fore,  a  Bi-static  operation  mode  could  be  considered 
and  also  multiple  locations  along  the  wall  can  be  consid¬ 
ered  to  perform  some  kind  of  spatial  filtering  to  reduce 
the  wall  clutter. 

VII.  Abbreviations 

DFT:  discrete  Fourier  transform;  SVD:  singular  value  decomposition;  STFT: 
short  time  Fourier  transform;  UWB:  Ultra  Wide  Band. 
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